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By L. M. K. Boelter, L. M. Grossman, 

R. C. Martinelli , and E. H. Morrin 


SUMMARY 


A critical comparison and summary is given of the various methods 
proposed to date for calculating the unit thermal conductance on the 
outer surface of a heated wing involving both laminar and turbulent 
boundary layers, and a new equation is suggested which should indicate 
the effect of the pressure gradient on the laminar heat transfer to a 
greater degree than do the expressions presented heretofore. For 
purposes of comparison the different equations are applied to a Joukowski 
profile for which the necessary data are accurately known and the results 
are plotted graphically. The unit thermal conductances in the laminar 
and turbulent regimes computed by the different methods are found to be 
in good agreement. A procedure whereby the equations for heat transfer 
from airfoil surfaces may be applied to a propeller shape is presented 
by means of an illustrative example. 


INTRODUCTION 


In the design of heating systems for wings for the purpose of thermal 
ice prevention a knowledge of the unit thermal conductances on the inner 
and outer surfaces of the wing is necessary. The internal conductance 
Vn=m probably a somewhat greater effeot on the value of the over-all heat 
transfer but the value of the external conductance is presumably the 
controlling factor in the case of the distribution of temperature along 
the wing. The external conductance is a function of many variables 
depending upon the form of the thermal and fluid boundary layers existing 
at the airfoil surface. ThiB report will be concerned solely with methods 
of c alcu lating this latter unit thermal conductance. Excellent summaries 
of the general problems involved in the design of wing heating systems 
are given in references 1 to 6. 
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Five general procedures^" have been suggested in the literature for 
computing the heat transfer into laminar boundary layers for incompressible 
flow along airfoil surfaces and four methods have been proposed for the 
turbulent regime. These are discussed in turn. In conclusion a new 
method, which is somewhat more complex than those previously published, is 
presented for the laminar regime and is intended to represent more nearly 
the aerodynamic and thermal conditions along a wing section. All the 
methods discussed treat the laminar and turbulent boundary layers separately, 
and to date no accurate means of locating the transition point has been 
evolved. Instead of an accurate specification of this .point, current 
practice favors the assumption that transition occurs at, or near, the 
point of minimum pressure. This criterion is of course inapplicable for 
the case of "laminar wings." 

The methods of calculation described in this report are each based 
upon one of four general principles (or simple modifications thereof). 

These are: 

For the laminar regime: 

(1) "Reynolds analogy" which states the equivalence of the 
equations for momentum and heat transfer at Prandtl modulus equal to 1 
and in the absence of a pressure gradient. 

(2) Pohlhausen*s exact solution of the differential equations 
for heat transfer in a laminar boundary layer along a flat plate for 
any value of the Prandtl modulus. Also, Colburn* s equation for heat 
transfer along a flat plate based upon Pohlhausen*s exact solution. 

(. 3 ) The postulate that the temperature and velocity distributions 
In the boundary layer are proportional, the factor of proportionality 
being calculated by means of a heat balance on the boundary layer 
(Squire* s method). 

(h) An integral heat balance over the laminar boundary layer. 

For the turbulent regime: 

(1) Reynolds analogy 

(2) Karman’s modification of Reynolds analogy 

( 3 ) Colburn’s equation for turbulent heat transfer along a 
flat plate which is Identical with the equations obtained from 
Reynolds analogy, except that the effect of the variation in the 
Prandtl modulus from unity is approximately accounted for 


^The paper presented by A. N . Tifford . (A.S.M.E. Aviation Meeting, 
Los Angeles, Juiie 19^*0 was not received before this report was written. 
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(It) An Integral heat balance over the turbulent boundary layer 

The methods .of analysis summarized in this report also illustrate 
several approx imat e procedures for. Indicating the effect of the 
pressure gradient existing over an airfoil surface upon the heat 
tr ans fer into the boundary layer. For purposes of comparison these 
procedures may be listed under four main headings: . „ 

(A) Substitution of the velocity near the airfoil surface 
(calculated by means of the pressure distributions about 
the airfoil ; see appendix B) into flat— plate equations and 
the substitution of approx im ate magnitudes of the local drag 

along the airfoil surface 

into heat— transfer equations derived from a consideration of 
flow over a flat plate 

(B) Substitution of the velocity near the airfoil .surface 
(calculated by means of the pressure distribution about the 
airfoil, see appendix B) into Colburn * b empirical heat- 
transfer equations for a flat plate 

(C) Heat balance on the boundary layer, including the effect of 
pressure gradient on the velocity distribution with the 
further postulate that the temperature distribution 1 b 
proportional to the velocity distribution 

(D) Heat balance on the boundary layer solving the heat transfer 
and hydrodynamic equations simultaneously 

With these classifications the methods described herein may 
be tabulated as follows! 
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The material in this report is divided into two sections; the first section describee the 
methods for calculations in the laminar regime and the second section describes the methods for the 
turbulent regime. 

The authors wish to record here their appreciation of the advice and assistance offered by 
Mr. G. Young during the preparation of this report. 

This work was conducted at the University of California under the sponsorship and with the 
financial assistance of ihe National Advisory Ccnmittee for Aeronautics. 
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A 

a 



J 

k 

Nu 

P 

lx 

E 

Ee c 


Ee D 


Ee, 


SYMBOIS 


area of airfoil equal to chord times spaa 


thermal diffusivity of fluid 



, square feet 
(sq ft)/(sec) 


drag coefficient 
lift coefficient 


( f D = |°fP Au cc 2 ) 
( f l = |c L pAu^ 


heat capacity of fluid at constant pressure, Btu/(lh)(°F) 

wing chord, feet 

diameter of equivalent cylinder, feet 

average unit thermal convective conductance, for length L, 
"between airfoil surface and air, Btu/(hr)(sq ft ) (°F) 

local or point unit thermal convective conductance between air- 
foil surface and air at any point x, Btu/(hr)(sq ftH^) 

local or point unit thermal convective conductance between air- 
foil surface and air at any angle cp, Btu/(hr)(sq ft)(°F) 

drag force, pounds 

lift force, pounds 


mechanical equivalent of heat (778), (ft— lb)/Btu 
thermal conductivity of fluid, Btu/(hr)(sq ft)(°F/ft) 


Nus selt modulus 
pressure, (lb)/(sq ft) 

heat transferred at any point x, Btu/(hr) 
ideal gas constant, (ft— lb) /(lb) (°E) 

Beynolds modulus based on wing chord 

Eeynolds modulus based on cylinder diameter ( ^ sa ^- 
Eeynolds modulus based on length x 
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temperature of fluid in boundary layer , °R 
temperature of fluid at edge of boundary layer, °E 
temperature of fluid in free stream, °R 
temperature of Surface, °E 

velocity of fluid near airfoil surface in x-direction calculated 
from Bernoulli’s equation, = — p (ft)/(sec) 

velocity of fluid in boundary layer in x-directicn, '(ft)/(sec) 

velocity of fluid in free stream in. x-direction, (ft)/(sec) 

velocity of fluid in boundary layer in y-direction, (ft)/(sec) 

length along airfoil profile measured from point of stagnation, 
feet 

distance normal to Airfoil surface, feet 
angle of attack of airfoil, degree 
dimensionless parameter in Pohlhausen solution 


thickness of laminar boundary layer, feet 
thickness of thermal, boundary layer, feet 
displacement thickness of hydrodynamic boundary layer 


Vo( 1 "») a 7 


feet 


displacement thickness of thermal boundary layer 

(i 4 *** 

"characteristic length" in turbulent boundary layer, feet 


angle between radius through, point on cylinder and radius 
through point of stagnation measured at axis of cylinder, 
degree; also functions defined by equations (^0) and (47) 
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P 
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7 
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To 


Er 


absolute viscosity of fluid, (lb)(sec)/(sq ft) 
kinematic viscosity of fluid, (sq ft)/(sec) 
weight density of fluid, (lb)/(ft3) 
mass density of fluid (?/g), (lb)(sec 2 )/(ftM 


momentum thickness of boundary layer 


a dimensionless quantity 





"eddy viscosity" defined by the equation 
(lb)(sec)/(sq ft) 

drag at the surface, (lb)/(sq ft) 


T _ _ du 

~ G ~ i 

P 3y 


Prandtl modulus 


^ 3600 uCpg ^ 


feet 


DISCUSSION OF METHODS 
Laminar Regime 


The first part of this sectionjconsists of a generalized discussion 
of the theory for the calculation of heat flow into 1 e-miner boundary layers 
The details of five methods are then stated at the end of the section. 

(jjEgynelds analogy .— With the usual postulates of boundary-: layer 
theory (see appendix B)7 the hydrodynamic equation for steady two- 
dimensional flow in a laminar boundary layer in the absence of a pressure 
gradient 2 is (reference 7, vol. I, ch. IV) 


3u 3u „ 3 2 u. 

3^ + T SP 


( 1 ) 


The equation for the temperature distribution in such a boundary layer is 
(reference 7, vol. II, p. 610) 

A flat plate oriented in the direction in which the fluid is flowing 
will -satisfy the requirements of equation (l) . ( pee appendix B.) 
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u 


ST 

3x 


+ V 


St _ v s 2 ! 
5y " Sr 


( 2 ) 


When the velocity and temperature of equations (l) and (2) are expressed 
in dimensionless fom by dividing, respectively, by u« the velocity of 

the free stream (in the oase of a flat plate the velocity in the free 
stream is equal to the velocity at the edge of the boundary layer) n-nri 
by AT = T 0 — Too, the difference in temperature between the surface and 

the free stream, there is obtained 


u 




v 



(3) 


and 


<£) >© 

+ V = -* 5— 

Sx Sy Pr Sy^ 


(M 


Inspection of equations ( 3 ) and (4) reveals that if Pr = 1 and if both 
equations have the same boundary conditions, the solutions of the two 

equations are Identical ^that is, — = and thus the temperature and 

velocity distribvitions are exactly similar. 


If the solutions are identical, then 



(5) 


or 



( 6 ) 


This fundamental relation expresses the well— known Keynolds analogy 
which states that there exists a direct relation between the temperature 
and velocity gradients, which respectively control the rate of heat transfer 
and the drag (shear) at each point on the surface of a flat plate losing 
heat to a stream of fluid. 
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Since 


Tq - H 


°f VP* 


( 7 ) 


y =0 


and 


To-^ 


7=0 


(7a) 


equation ( 6 ) may be written as 


TCtjII, 


Pr 


Oftr 


( 8 ) 


E “ 


Since in the derivation of equation ( 8 ), Pr was postulated to he equal 
to unity, the equation should he written as 


f 03C = £fj 
2 


(9) 


It cannot he overemphasized that equations ( 6 ) and ( 9 ) are valid 
only if each of the two fundamental conditions upon which they Eire based 

is satisfied, namely that ^ = 0, and Pr = 1. If the pressure gradient 

along the surface Is not zero, equation (1) no longer describes the 
velocity In the bo und a r y layer hut a term Involving the pressure gradient 
must he added. That is, for laminar boundary— layer flow over a surface 
other than a flat plate, the hydro dynamio equation is 


u 


Su 

dx 


+ V 


du _ 1 dp 

§y p 


+ v 


3^u 

dy 2 


(10) 


Instead of equation (l), and the dimens ionless equations for the velocity 
and temperature will no longer he identical. Reference to equations ( 3 ) 
and (4) indicates clehrly also that unless Pr = 1, the equations will no 
longer he identical. 

A simultaneous solution of equations (10) and (2) is difficult except 
for certain special variations of with x (reference 7 , vol. II, 

pp. 63 I— 635 ) (appendix B). These special cases indioate, however, that an 
approximate correction for the pressure gradient existing along the 
airfoil surface may he made by 



10 


MCA TN No. 1453 


(a) Substituting U, the velocity near the airfoil surface 
(calculated freer the pressure distribution about the airfoil) 
for u M in the Reynolds analogy , 

(b) Substituting the local drag coefficient along the airfoil 

(based on ^ along the airfoil) for the flat-plate 

drag coefficient required by the Reynolds analogy 

Thus the Reynolds analogy modified approximately to account for pressure 
gradient is 



(H) 


( 2 ) Pohlhausen solution .— Pohlhausen (reference 8) has solved equations (l 
and (2) simultaneously for magnitudes of Pr other than unity by substituting 
the Blasius series solution for the velocity in equation (2) to obtain the 
temperature distribution in a laminar boundary layer along a flat plate. 

His solution may be written 

IBS * 

H=o 


0.332 ^PF (T - T ) 

V vx 


( 12 ) 


Since the Blasius solution for the velocity distribution in a laminar 
boundary layer along a flat plate is (reference 7, vol. I, p. 135 ) 



0.332 



u_ 


( 13 ) 


the ratio of the temperature gradient at the wall to the velocity gradient 
now becomes 

Wy-o = pj.l/3 do - M 

(jg) ^ 

\Oy/y=0 

When this expression is compared with equation (6) 

Pohlhausen solution yields the same results as the 
except for the inclusion of the term involving the 
equation (14) may be utilized to calculate heat losses from a flat plate 
for magnitudes of the Prandtl modulus differing from unity. The condition 
that the pressure gradient along the plate is zero must still be fulfilled 
however, and equation (14) is again strictly valid only for a flat plate. 


( 14 ) 


it is noted that the 
Reynolds analogy 
Prandtl modulus. Thus, 
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As before, equation ( 14 ) may be rearranged, for 


T o 



G f v P u « 2 

2 


and 





Thus by substitution into equation ( 14 ) 



360070^ 




(15) 


Colburn has demonstrated that equation (15) satisfactorily correlates 
expe ri mental data for the heat loss from a flat plate. 

By the use of the equation for the local drag coefficient along a 
flat plate (reference 7, vol. I, p. I35) 


Of, - 0.332Ke x -^ 


(16) 


Equation (15) may be rewritten as 

Pr 2 /3 




360070^ 


0.332Re x 


-0.5 


(IT) 


where 


Ee x 


US 


Equation (15) may be modified as before for application to airfoils by 
the substitution of U for u w and the calculation of the local drag 

coefficient for the airfoil. Thus equation (15) bee ernes 




36007C p U 



(IS) 
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Equation (i.7) may also be modified for application to airfoils "fay the 
substitution of U for u*, . The correction for, differences in the local 
drag between the airfoil and the flat plate is not accounted for, however. 
Thus 


V.- 

36007C p U 


0.332R6X -0,5 


(19) 


Equation (19) is termed the Colburn laminar e quia t ion in this report. 

( 3 ) Proportionality between velocity and temperature distributions .— 
If it is postulated that the temperature distribution is everywhere 
proportional to the velocity distribution, then 

= 1 _ It ' 


If, in addition, the ratio — is expressed as a known function of the 
V 

parameter — where is the hydrodynamic-boundary— layer thickness, 

S 1 

that is. 



then from, the similarity of the temperature and velocity profiles 


T 

T 


- T„ 


- T 0 


= 1 - 



where f is the same function as in the hydrodynamic case and Br, is the 


thickness of the thermal boundary layer. The problem of determining the 
temperature distribution then reduces to one of expressing $2 in terms 

of the known value of Si. The ratio is in general a function of both 

6 1 


the Prandtl modulus and the pressure distribution and reduces to unity 
only in the case of heat transfer free, a flat plate at Pr = 1. Squire's 
contribution consists in deriving this functional relation by means of the 
energy balance described in the following section. 


(4) Heat balance .— When an incompressible fluid _is considered and the 
'’dissipation" effect is neglected, an energy balance may be made across a* 
section of the boundary layer of width dx at a point x, (see fig. 1) 
as follows. The amount of excess heat being carried across a nor mal 

0 00 

to the surface at x is Cp? I u(T — T«J dy and the difference 
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between the value of this expression at x + dx and the value at x 
must he equal to the amount of heat entering the fluid from the surface. 


that is, — k( 


/ y=o 


For uniform density and zero dissipation, the energy 


equation for the thermal layer is 


& 


u(T - ^o) dy 



(18a) 


If u and T are known, the point unit conductance along any 
surface as a function of y is then given by 


f Qx _ i a 

36007CpN V Sx 


U 



( 19 a) 


It should be noted that, in general, principles (l) and (2) are based 
on a knowledge of the velocity gradient at the surface of the body, -whereas 
principle (3) requires that the velocity gradient at the wall as well as 
the velocity distribution be known. Principle ( 4 ) requires the knowledge 
of both the temperature and- velocity distributions in the laminar boundary- 
layer but does not require an exact knowledge of the velocity and temperature 
gradients at the wall. 

In the following paragraphs several methods which have been proposed 
for the calculation of laminar-boundary— layer unit conductance along wings 
are compared, 

Method of Allen and Look .— Allen and Look (reference 9 ) treat only 
the case of heat transfer into laminar boundary layers, basing the 
derivation upon Reynolds analogy as expressed by equation (ll), 

fat Est 

36007CpU 2 


or 


/ 

f c x - |c fx 7C p U X 3600 


In order to evaluate the local drag coefficient C-p^, the Blasius 

■type velocity profile, which represents the solution of the differential 
equation for flow in an incompressible boundary layer along a flat plate 
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(with zero pressure gradient), is postulated by the authors to apply over 
the entire airfoil. In terns of the laml nar-boundary— layer 3- thickness 8, 
the local drag coefficient along a flat plate hoc canes 

(») 

When equation (20) is substituted in equation (11) and rearranged, 
the following equation is obtained 

5^ = 0.765 = 0.765Pr 


or 


f _ 0*765 vp,, _ 0.76f5k 

c x “ 8 6 


( 21 ) 


beoause the Prandtl modulus was assumed to be unity . in order to obtain 
equation (11). 

Equation (21) suggests that heat is conducted through the boundary 
layer by conduction across the thickness 5. Physically this is not 
the case; heat is transferred into the boundary layer and carried within 
the boundary layer along the surface of the plate. Ho heat is transferred 
to the free air stream, by transfer across the boundary layer. 

In order to account approximately for the effect of the pressure 
gradient existing about the airfoil, the thickness of the laminar boundary 
layer 5 Is computed by substituting the Blasius distribution in the 
von Earmin momentum condition and integrating. This operation yields an 
equation for 8 2 of the form (reference 10) 



Equation (22) for the boundary— layer thickness represents the point 
in the derivation of Allen and Look where the treatment differs fraa the 
case of heat transfer in a laminar boundary layer along a flat plate. 

'^The thickness of the laminar boundary layer 8 is arbitrarily defined 
as the distance from the solid-fluid interface to a point in the boundary 
layer where the dynamic pressure is one— half of its value outside the 
boundary layer. (See appendix B) 
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of flow along a flat plate ( bhat is. 


has been used, substitution 


Although the Blasius velocity profile which is only applicable to conditions 

s = °) 

of the equation for the Blasius velocity profile into the m ome ntum, equation 
containing a term involving the pressure gradient should yield a value of 5 
more nearly representative of conditions at the surface of an airfoil. 


0 


Combination of equations (21) and (22) yields the final equation: 


f cx . 0.332 
3600C p 7lT 



(23) 


The equation is written in this manner to allow ready comparison with 
flat— plate heat-transfer equations. The Prandtl modulus does not appear 
in equation ( 23 ) because it was initially postulated to be unity. 

Method of Frick and McCullough .— The method of Prick and McCullough 
(reference 11) for treating the heat transfer through laminar boundary 
layers is similar to that of Allen and Look but is generalized to include 
values of Pr other than unity by utilizing Pohlhausen’s exact solution 
of the differential equations for heat transfer in a laminar boundary 
layer along a flat plate. 

As discussed in the introduction, the Pohlhaixsen solution yields the 
expression: 




As in the method of Allen and Look the velocity gradient at the wall is 
de lived from the Blasius solution for a flat plate for the condition that 

at y = 6, = O. 707 . (See footnote 3«) 


Thus 


therefore 




= O .765 




(24) 


(25) 
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Thus 

' X s " K^y-O = fc *f T ° " = °* 765 ^ l (T 0 " T “) (26) 

or, for a magnitude 14 ' of Pr = 0.760 

?c x = 0-700 | ( 27 ) 


where 5, as in the Allen and Look report, is given "by equation (22). 
The substitution of this value of 5 then partly accounts for the flow 
conditions along surfaces other than flat ones. 


The difference between the 'coefficient 0.700 as found by Frick and 
McCullough and the value of O .765 found by Allen and Look arises from, 
the presence of the oube root of Pr in equation ( 26 ). 


In general, for other magnitudes of Pr, by substitution of 
equation (22) into equation (26) 


f °x 

3600C p 7U 


0.-212 



J 

M ^ 

c~ O 

£) 8 - 17 <!) 

1 

Hi) 


- 1/2 


( 28 ) 


Basically this solution is identical to that of Allen and Look, with the 
exception of the use of the Pohlhausen solution instead of the Keynolds 
solution for heat transfer from a heated plate. 

Method o f Martlnelli, Boelte r, and others .— Probably the simplest 
approximate method of computing the point unit thermal conductance of 
the laminar boundary layer over wings is that of Martlnelli, Guibert, 
Morrin, and Boelter (reference 13). The airfoil surface is conceived 
as a combination of a cylinder close to the leading edge and a flat plate 
beyond; the known data and equations for the heat transfer from these 
surfaces are then applied to the ideal surfaces, 

(a) Near stagnation point 

For the magnitude of the unit thermal conductance at the 
stagnation point of a cylinder. Squire's analytical solution of the 

s ummar y of data (reference 12 ) reveals that the value of Pr 
for the temperatures usually encountered in wing anti— icing is closer to 
0.72. A magnitude of Pr = O .72 is thus utilised in all laminar airfoil 
calculations in the present report. 
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differential equations for the boundary layer (reference 7, vol. II, 
p. 63 I) yields 


Nu . = 1 . l4Pr 0 • ^e-n 0 ’ ^ 

stag v 


(29) 


The point conductance along the leading edge at any angle (up to about 
9 = 9 O 0 ) measured frcm the stagnation point (reference 14) can then be 
approximated by the equation 


f 



1 . l4Pr ^ ‘ • 5 


1 



(30) 


or if the properties of air are expressed by a power function in the 
absolute temperature T 



(b) Remainder of wing (laminar flow) 

For the heat transfer along the portion of the airfoil section 
which is considered to behave as a flat plate, the equation of Colburn 
based upon Pohlhausen’s analytical solution is used, (reference 15 ) 


^c- 


36000^711 


(Pr ) 2/3 = 0.332ESX " 0 ' 5 


(19) 


where U Is the velocity near the airfoil surface at the point x, 
computed frcrn the pressure distribution existing about the airfoil. 
(See appendix B.) The coordinate x is measured from the stagnation 
point. When the properties of air are expressed as a function of 
temperature 


f Cx = 0.0562 t0.50(3£)°-° ( 3 2) 

For comparison with the previous two methods, equation (19) nay be 
written as 


XT* 



36 ooc P 7 U 


. (pr r 2/3 


(33) 
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The heat— transfer equation (3.0) used in Martinelli ’s method for 
calculating the point conductance near the forward stagnation point is 
obtained from analysis and data on heat transfer near the stagnation 
point of a cylinder. The equations used beyond the leading edge in the 
method are, of course, purely flat— plate relations, but the substitution 
of the actual velocity U at the edge of the boundary layer for u^ is 
intended to account approximately for the effect of the pressure 
distribution along an airfoil surface. 

Method of Squire .— The procedure adopted by Squire (reference 16) 
for the calculation of the point unit thermal conductance over the outer 
surface of an airfoil involves more lengthy computation than any of the 
preceding methods, but probably represents the most rational analysis of 
the aerodynamic and thermal relations thus far published. 

It can be shown that when the Peclet modulus (ticd/a) is large, the 
conditions with respect to the temperature distribution near the surface of 
a heated body past which fluid is flowing are the same ae those with 
respect to the velocity distribution when the Reynolds modulus 1 b large. 
That is, a "thermal boundary layer" of small thickness exists near the 
surface, in which the temperature falls rapidly frcan its value at the 
surface to the temperature in the free stream. Thus it follows to define 
a "thermal displacement thickness" 6 2 for the thermal boundary layer 
in a manner analogous to the hydrodynamic case, (appendix B) by the 
equation. 


5 


2 



( 3*0 


in which T, T 0 , and are the temperatures in the bou nd ary layer, at 

the surface, and in the free stream, respectively. 

In Squire's derivation an hypothesis of fundamental importance is 
made; namely that the temperature distribution is of a similar form, to 
that of velocity. It should be clearly borne in mind that when this 
postulate is made it is assumed that the pressure gradient affeots the 
temperature distribution in exactly the same manner as it does the velocity 
distribution. This is not strictly correct, and the accuracy of the 
approximation is still open to experimental verification. (See appendix C.) 

With this hypothesis, utilizing the Blasius series solution for the 
velocity distribution in the laminar boundary layer of a flat plate, there 
is obtained 




( 35 ) 
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or 



5 i w 


06 ) 


The equation for the hydrodynamic—' boundary— layer displacement 
thickness (appendix B) may he written as (reference 17 ) 



2.960 

T-vx 



(37) 


so that 


or 




3600C p tU 



(38) 


(39) 


The ratio 


8 1 , , 

■== appearing in equation ( 39 ) is a function of U and. 
5 g 


Pr and was obtained by Squire by means of the integral heat balance 
derived in section entitled "Heat balance." When the Blasius distribution 
for velocity a n d temperature is substituted in this heat balance & 2 W 

be expressed in terms of the known value of 5^ by the relation 



( 40 ) 
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where <p is a known function of its argument which Squire tabulates for 
S P 

values of — between 0.5 and 2-.0. If the pressure gradients are not too 


83 

large, the ratio — - very nearly equals Pr' 


.1/3 


It is important that one recognize wherein Squire's method is similar 
to the approach used by Look and Allen, Frick and McCullough, and Martinelli 
and in what way it differs from, these. It will be recalled that in the 
latter three methods the authors use the hydrodynamic equations and 
heat— transfer equations for a flat plate and include the effect of the 
pressure gradient only in computing the value of the boundary— layer thickness, 
which appears in their final equations for the local unit thermal conductance, 
and by substituting the actual velocity at the edge of the boundary layer. 
Although it is true that Squire also used the flow and heat— transfer 
equations which are valid only for a flat plate, his method represents 
an advance in that it uses a corrected value of the thermal— boundary— layer 
thickness computed as a function of the hydrodynamic-bcundary— layer 
thickness in the final equation for f c . 


Method of this report .— It will be noted that each^ of the preceding 
methods involves a camp remise or combination of hydrodynamic and heat- 
transfer equations which are strictly applicable only in the case of 
boundary layers along plane surfaces, together with the momentum equation 
which is valid for any curved surface up to the point of separation, to 
yield final expressions for the unit thermal condii.ctance over an airfoil 
shape. 

The authors of the present report have attempted an analysis of the 
problem based on relations true for any boundary layer regardless of the 
shape of the section and without recourse - to the Blasius flat— plate 
solution. 


Two such general expressions are available: one being merely an 

integrated heat balance across the boundary layer, andthe other an 
integrated force balance. If one postulates -steady flow, an incompressible 
fluid, and no dissipation of kinetic energy, the fir sty often called the 
"energy equation" of the boundary layer, is 


Bx 



uT dy — T^ 


Sx 



u dy = — 



(>» 1 ) 


and the second, is the well.— known von Kazman momentum equation (appendix B) 


u dy' 


P Sx _ Vdy/^o 


(■’> 2 ) 


ox 


'6 


0 


u 2 dy - U 

dx 


0 
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When a fourth power polynomial Is postulated to represent both 
velocity and temperature distribution in the boundary layer, the 
following equation is obtained with the use of the proper bou ndar y 
conditions. (See appendix A.) 


u 

U 



where 



Equation (4l) may be written in the following foim: 




, 5 Hu 




, 5 Eu 


JO - XI - ox Jo U 


= dy = -a 


(g). 


y=o 


which, upon substitution of the polynomials for the velocity and 
temperature distribution and integrating, yield n 


^ 


where 


and 


<P = [(-0.134 - 0. 01H)Pr —1 /3 + (0.021 - 0 . 006 X)Pr “ 1 

+ (- 0.006 + O.OOllJPr”^/ 3 ] 


5h ~ -n -1/3 , 

“g ~ (See appendix A. ) 


the 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 



22 


HACA TR Ho. 1453 


Therefore 


f Cx - 3600rC p ^ (USg?) 


(49) 


Details of the method of derivation and the procedure for calculating 
the point conductance in this method will he found in appendix A. 

It was previously pointed out that Squire’s method represented a 
oloser correspondence to physical, facts than those previously proposed 
in that a general energy balance on the boundary layer was used to 
calculate the ratio of the thermal— to the hydrodynamic— boundary— layer 
thickness. The method does, however, employ the Blasius series solution 
for the flat plate to represent the velocity distribution in the boundary 
layer. In the aforementioned method this latter approximation is 
eliminated by using both fluid-flow and heat— transfer equations which are 
general for the laminar boundary layer along any type of section and which 
take into account the pressure distribution actually existing over the 
airfoil surface. Specifically, the Fohlhausen polynomial 1 b used in 
preference to the Blasius series to represent the velocity distribution 
in the boundary layer. 

In so doing, the limitations of the Fohlhausen method should be 
clearly borne in mind (appendix B) . It is generally conceded that this 
method gives a good representation of the velocity distribution in a 
ldminar boundary layer in which the fluid is being accelerated and is 
fairly accurate in regions of retarded flow at positions distant from 
the point of separation. The method bee ernes less and less accurate as 
the point of separation is approached and breaks down completely upon 
reaching that point. 

Summary of methods .— The final equations for the laminar point 
unit thermal conductance for each of the four methods previously 
published are as follows : 


Allen and Look: 



Frick and McCullough: 

f °x = 0.332 py-2/3 
3600C p 7U \fnex_ 



(23) 


(28) 
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Martinelli and others: 


f °s = 0.332 ^-2/3 
3600C p 7U \fR e^ 


Squire: 



3600C p 7lJ 



( 33 ) 


(39) 


Comparison of equations (23), (28), (33), and (39) reveals that: 

(1) Equations (28), (33), and ( 39 ) account for variations of Pr 
from unity, hut equation ( 23 ) does not. 

(2) All four equations account approximately for the pressure gradient 
existing along the airfoil hy substituting the velocity at the edge of 
laminar boundary layer for Uoo in flat— plate relations but, in addition, 
equations ( 23 ), ( 28 ), and ( 39 ) make further approximate corrections by 
means of the terms in brackets. 


( 3 ) Equations (28) and ( 33 ) are identical for the case of heat transfer 
over plane surfaces, and all the equations are identical for heat transfer 
from a flat plate at Pr = 1. 


The next section of this report is devoted to a discussion of heat 
flow into turbulent boundary layers. A generalized discussion of theory 
is followed by details of four calculation procedures. 


Turbulent Regime 

Explorations of the velocity immediately adjacent to streamlined 
surfaces along which a turbulent boundary layer exists are, not numerous, 
but measurements of the vplocity near the walls of tubes in which a fluid 
is flowing turbulently have been accurately performed (references 18 to 20 ) 
and reveal that a completely turbulent boundary layer probably does not 
exist. Near the walls of the tube there is found to be a laminar "sublayer" 
in which the flow remains vise crus even though the fluid far frcm the surface 
flows turbulently. For purposes of analysis a transition layer, sometimes 
called the "buffer layer" may be visualized as existing between this 
laminar sublayer and the turbulent fluid. In the laminar sublayer viscous 
forces predominate ; in the turbulent region "eddy" forces are controlling; 
whereas in the buffer layer both viscous and eddy forces are of the same 
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order of magnitude. Reynolds (reference 21) postulated that in the 
turbulent region, momentum and heat are transferred hy similar mechanisms. 
In the laminar sublayer, molecular heat transfer occurs, whereas in the 
buffer layer both moleoular and "eddy" heat transfer take place. 

On the basiB of these concepts, a boundary— layer heat-transfer 
equation and a boundary-layer hydrodynamic equation have been written 
(reference 7 > vol. II, p. 650). The heat— transfer equation is 


u dT, v dT_d/v 3 t 
S y §y( Pr §y 


+ e H 


ST\ 

sy ) 


where gjj is the so-called "eddy diffusivity" for heat. 
The hydrodynamic equation is 


u 


5u 

S3E 



+ e 


Su 


) 


(50) 


(51) 


where is the so-called eddy diffusivity for momentum. The velocities 

u and v in equations (50) and (51) are mean values with respect to time. 


On the basis of the momentum-transfer theory exact- similarity between 
the temperature and velocity distributions therefore exists if 


(1) The eddy diffusivity for heat eg equals the eddy diffusivity 
for momentum ey. (This statement is the original basis for the 
Reynolds analogy as developed by Stanton (reference 22).) 


(2) The Prandtl modulus equals unity. (This statement Is necessary 
beoause of the existence of the laminar sublayer and buffer layer, 
in which viscous stresses and molecular heat transfer are important.) 


(3) The pressure gradient 


ctx 


is zero. 


If all these conditions are satisfied and if the temperature and velocity 
distributions have the same boundary conditions,' the solutions for velocity 
and temperature will be identical. Thus, for turbulent flow, as in the 
case of laminar flow: 



( 52 ) 
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As in the c&se of a laminar boundary layer, this equation can he rewritten 
as: 


or since Pr = 1 


3600C-.7UCO 



360CX3pru 0o “ 2 


(53) 


. (5M 


This, equation is usually referred to as Reynolds analogy and , as 
■written, applies only to heat transfer from plates = 0 ^ at a uniform 

temperature to a fluid with Pr = 1. Modifications of the analogy to 
apply for other magnitudes of Pr and for flow conditions in which a 
pressure gradient exists will he discussed helow. 

Modification of Reynolds analogy to account for pressure gradient 

the pressure gradient is not zero, then equations ( 50 ) and ( 51 ) are no 

longer analogous and the exact equivalence of u and T, even for Pr = 1, 
breaks down. 


Calculations for an airfoil on the basis of a heat balance in which 
the temperature and velocity distributions in the turbulent boundary layer 
are assumed to vary identically with pressure gradient (reference 1 6 ) 
reveal that the pressure gradient may have a large effect on the rate of 
heat transfer. The calculations presented in reference 1 6 , however, ' 

probably overemphasize the effect of ^ because actually, the pressure 
gradient influences the velocity much more than it does the temperature, 
as noted by the fact that the term 


enters the velocity equation 


directly, but influences the temperature only through its effect on the 
velocities u and v. 


Thus the exact influence of on the rate of heat transfer is 

ox 

difficult to establish; one may probably say that the effect should not 
be pronounced, but rather of secondary importance. For lack of a 
simultaneous solution of the heat— transfer and hydrody nami c Equations 
for the turbulent boundary layer, including the effect of pressure gradient, 
approximate corrections to the Reynolds analogy (similar to those made 
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for the laminar boundary layer) to account for ^ have been presented 

in the literature* In general., these approximate corrections consist of 
substituting U, the velocity near the airfoil surface (calculated by 
means of the pressure distribution existing about the airfoilj see 
appendix B), for u*, in equation ( 54 ) and calculating the drag coefficient 
Cf x for the airfoil at the point in question. The Bey no Ids analogy then 

becomes: 


36coc p ?u 



( 55 ) 


Modification of Reynolds analogy for magnitudes of Pr other than 
unity.— Since the Prandtl modulus for air is less than unity, the Eeynolds 
analogy must be modified to allow calculations for other magnitudes of 
the Prandtl modulus. Several methods are available, all being based on 
a n alyses of heat transfer to fluids in tubes: 

(a) Yon Kannan modification: By analyzing the heat transfer from, a 

tube to an enclosed fluid flowing turb ulently , von Kazman obtained 
(reference 23): 


C f 



(b) Boelter, Martins ill, and Jonassen modification: By extending 

the von Karrnan analysis to include a more precise consideration 
of the turbulent region, Boelter and others obtained for flow 
in tubes only (reference 24 ): 


Cf 


Zffi. 


max 

'mean 


3600Q v 7V m 


( 57 ) 


Pr + logg (1 + 5Pr) + 0.5 logg 


Ee 


60 \| 


Cf 


(c) Colburn, by empirical, correlation of data on heat transfer from 
flat plates, obtained (reference 15): 


3600C p 7U m 




( 56 ) 
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numerical calculation of the von Karman and the Boelter method 
reveals that for a range of Pr "between 0.5 and 10 , the complex expressions 
involving Pr reduce to that of Colburn with fair accuracy (reference 24). 
Thus, for fluids with 0.5 < Pr <10 /the modified Reynolds analogy may "be 
rewritten as: 


3600C p7 U 


Pr' 


2/3 = 

2 


(59) 


When the equation for the local drag coefficient along a flat plate 
(reference 7, vol. II, p. 362 ) 

=, o. 0296 Ee x “°' 2 ( 60 ) 

is substituted into equation ( 59 )* the expression becames: 

3^7? - O.oasfipr-^/^- 2 (61) 

This equation Is called "Colburn’s" turbulent heat—- transfer equation 
In the remainder of this report. 

Heat balance .— If the velocity and temperature distributions are 
accurately known In the boundary layer, a heat balance will yield a 
simple method of obtaining the unit thermal conductance frcm the airfoil. 
As in the case of laminar flow 


fc x 

3600Cp7U 



U Sx 



(T - Tj 
(*o “ 



(19a) 


The difficulty in application of this method lies in the necessity 
for accurate knowledge of the velocity and temperature distributions. 
The methods based on the Reynolds analogy and its modifications require 
only a knowledge of the velocity gradient at the surface of the solid, 
as may be seen frcm. an inspection of equation ( 52 ). 

Method of Frick and McCullough .— Utilizing the Reynolds ana logy. 
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Frick and McCullough, obtain 


SE - n t 

T _a P c 


ph*) 


(63) 


In order to appiy this equation to heat transfer from airfoil surfaces, 
the local shear r Q is calculated by means of the von Karman expression for 

the skin friction over a flat plate with a fully developed turbulent 
boundary layer. Thus 


To - 


pU‘ 


(64) 


where 


£ = 2.557 log e ^4.075 


( 65 ) 


in which the momentum thickness 0 is calculated from airfoil boundary- 
layer data. (See appendix B.) Substituting for t q in equation ( 63 ) 

and rearranging yields 




(66) 


If one considers 


& 


Be, 


a characteristic length for the turbulent 


boundary layer, then 


CL X k 

X = ^ 


Pr (%o - T 0 ) 


(67) 


or 
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Letting ^ px - 0.760, Frick and McCullough obtain^ 



0.760k 

®T 


( 68 ) 


This characteristic length should not he confused with the thickness 

of the turbulent boundary layer. For a flat plate this latter quantity is 
given by the expression (reference 7 * vol. II, p. 362) 


- = 0 . 37 Ee x "°* 2 
x 


(69) 


from, which it is seen that &ji is approximately proportional to x and 
inversely proportional to the boundary— layer thickness. 

In summary it may be stated that the treatment of turbulent-boundary— 
layer heat transfer in the report of Frick and McCullough assumes complete 
equivalence of skin friction and heat transfer and uses the best method 
known for the calculation of turbulent skin friction. 

In order to compare the available methods for calculating heat 
transfer into turbulent boundary layers, equation (68) is rewritten 
for ready comparison with the flat— plate equations to be derived in the 
next section. The equation of Frick and. McCullough becomes : 

— = 0.0296Re x “°* 2 ( - 2^5 -) (70) 

3600Cp7C \0. 02965 / 

The Prandtl modulus does not appear in equation (70) because a magnitude of 
Pr = 1 was tacitly assumed in its derivation. 

As in the case of the laminar sublayer Pr =* 0.72 is more correct 
and is utilized in the calculation of this report. 

^It should be kept in mind that equation (68) is based on equation (62) 
Equation ( 62 ) does not include Pr, but in its derivation a Pr value of 
unity is postulated. In order to obtain equation (63) from equation ( 62 ) 
a Prandtl modulus of unity is utilized. In order to obtain equation (67) 
from equation (63), equation (63) was multiplied and divided by Pr Thus, 
as long as proper values of n, Cp, and k are utilized in equation (67), 
substitution of the proper magnitude of Pr, even though differing from 
unity, will yield a value of f 0 which in reality Is calculated for 

Pr = 1 and is therefore too lew. In a later section the equations of 
Frick and McCullough are modified to correct approximately for magnitudes 
of Pr not equal to unity. 
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Method of Martins 111 and others .— For the case of turbulent flow 
beyond the leading edge, the Colburn equation for turbulent heat transfer 
from flat plates is utilized. Thus 

— ^ = 0.0296Pr -2 / 3 Ee x -0 * 2 (71) 

3600C p rU 

or expressing the properties of air approximately by a power function 
in T the absolute temperature. 



0.5H 


P.3 



(72) 


When equations (70) and (71) are .compared two points of difference are 
noted: 

(1) The tem in parentheses in equation (70) represents an approximate 
correction to the flat— plate equation for the pressure gradient 
existing about the airfoil, in addition to the use of the velocity 
near the airfoil surface in the flat-plate Reynolds analogy. 

(2) Equation (70) does not involve Pr, whereas equation (71) includes 

Pr — 2 /^. Since the latter equation accounts for the Prandtl modulus 
more accurately than the former, the equation of Frick and 
McCullough, to be more correct, should be written as: 


* = 0.0296Pr~ 2 / 3 Ee x " 0,2 

3600C p7 U 


Be * 0 * 2 \ 
0.0296^/ 


(73) 


Magnitudes of f Gx calculated from the original Frick and McCullough 

equation are thus too low by the factor Pr“ 2 /3 # Equation (73) is called 
the "modified Frick and McCullough equation" in this report. 


Methods of Squire .— Squire (reference 16) presents two methods of 
calculating f c for an airfoil. The first method employs the Reynolds 
x / ✓ 

analogy as modified by von Kazman. Thus 


36000^,7X1 




2 


(7*0 


1 + 5 \| 2 \ (Pr " X) + lo8e 1 + 5 (Pr “ 1} 


°f„ 


is calculated at each point along the airfoil by utilizing the 


where 



von Barman expression for skin friotion along a flat plate with a fully developed turbulent 
boundary layer, as outlined in the method of Frick and McCullough, 

In order to compare equation (74) with those for flat plates. It may he rewritten as: 


—ISx 

3600C p7 U 


0.0296Rq x “°' 2 



-L 

|(Pr - 1) + logg £ 


Hex 


0.2 


1 + 6 {Er - 1] |/V' 


02965 ' 


(75) 


The first tern in parentheses corrects the Reynolds analogy for the flat plate for 
magnitudes of Pr other than unity. RumerioaUy (for small values of Pr) the term 1 b 

approximately equal to Pr* 2 ^. The last term in parentheses represents the approximate 
correction to the flat-plate equation for the pressure gradient existing about the airfoil 
and is identical to that utilized by Frick and McCullough. 


To indicate the probable effect of pressure gradients. Squire presents an alternative 
form of Reynolds analogy by assuming that the tenqperature and velocity distributions are 
exactly similar in the presence of such gradients. 

The energy equation for the boundary layer is 


A 

Sr 


u(T - T„) dy 

Jo 


3600C p 7 


(76) 


When it is assumed that the temperature and velocity distributions are similar 
(Pr = 1) 


T - T w = x _ u 
T - T U 

O °o 


(77) 


u> 

W 


i 
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equation ( 76 ) 136001168 


f C-v 


3600C p 7U 



de u 

se + tr 


t 


( 78 ) 


where 0 is the momentum thickness of the boundary layer and equals 


2 ( 1 - 2 ) 

0 M u/ 


dy, (see appendix B) and the primes denote differentiation 


with respect to x. From the momentum theorem, however. 


|£ t 3 l£ 5 _uLef|- + 2 N ) 

a* p u 2 u V s / 


(79) 


so that 


fc x 

3600C p 7U 


W 





( 80 ) 


As has, been mentioned previously, a plot of £ 0 ^ against x/c for the 
usual Reynolds analogy and for this latter modification indicating the effect 
of shows large deviations for an airfoil example at — > 0.4. 


Because 


W 


equation ( 80 ) then directly indicates the 


comparison between the additional correction for pressure gradient and 
the usual Reynolds analogy for a flat plate. 


Also since the momentum theorem (equation (79)) s-nd the energy 
balance are equally valid for laminar and turbulent boundary layers, 
equation ( 80 ) applies as well to the laminar case, and the term to the 
right of the minus sign affords an indication of the relative effect of 
pressure gradient on the heat transfer into laminar and turbulent 
boundary layers. 


In order to compare equation ( 80 ) with the turbulent equation of 
Frick and McCullough and the flat— plate relations, it may be written In 
the form 



fc. 


^(YL.-yTT 
^ P / ~ 


^ W * vu 7 u ^ j L \ / J 


( 81 ) 


summary of methods .-- The final equations for the point unit thermal conductance for 
turbulent boundary layers, are as follows: 


Frick and bfcCullough (method I) 


36ooc p7 U 

Modified Frick and McC ull ough (method U) 


^ /> I Ua 

0.029^^^' 


An/ T5rt 0.2 \ 

\^e) 


-pL- - 0 ( -fo - l 

36ooc p7 U U.0296C 2 , 


na VH nnrl n+tia™ • 


( 70 ) 


( 73 ) 
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3600c 7 U 
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Squire (method I) 
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36000^ 


0.0296^ 
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Squire (method. II) 



The bracketed, terns in these equations directly indicate the corrections 
made for pressure gradient and variations in the Prandtl modulus. 

Comparison of equations (70), ( 73 ) * ( 71 ) * ( 75 )* and (8l) reveals 

that: 


(1) Equations (73)* (71)* and (75) account for variations of Pr from, 

unity, but equations (70) and (8l) do not. ' ~ 

( 2 ) All five equations account approximately for the pressure gradient 

existing along the airfoil by substituting the velocity near the airfoil 
surface TJ for in the flat-plate relations. In addition, equations (70), 

(73), and (75) make further approximate corrections for the variation of 
point drag coefficient along the airfoil surface. Finally, equation (81) 
includes a further corrective tern which results from the calculation of 

a heat balance on the boundary layer. The last equation probably over- 
emphasizes the role of the pressure gradient. 

(3) Equations (70) and (71) are identical for heat transfer from a 
flat plate and have been checked experimentally for this case. Equation (75) 
is practically identical to equations (71) and (73) for heat transfer 

over a flat plate, and basically probably accounts for variations of Pr 
from unity more precisely than equations (71) and (73) which are based on 
an empirical correlation of experimental data. 

(4) Although equations (70) and ( 71 ) are identical when gjj -= 0 , the 

results from these equations are strictly applicable only to fluids with 
Pr - 1 . 

(5) All equations are identical for heat transfer from a flat plate 
to a fluid with Pr = 1. 


DISCUSSIOW OF WUl'QSHICAL EXAMPLES 
Methods Employed in Laminar Regime 


In order to compare the various methods described for calculating the 
laminar point unit thermal conductance over wings, an airfoil shape was 
selected whose aerodynamic characteristics are very accurately known. 
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Such a section* for which the pressure and velocity distribution has been 
calculated theoretically and thoroughly tested experimentally* is a 
Joukowski profile whose characteristics are shown in figure 2. This 
particular profile is called "aerofoil A" by Bairstow and numerous tables 
of data concerning it will be found in reference 25 . 

For purposes of calculation* an airfoil with a chord of 7-78 feet 
was chosen which was to be maintained at a constant temperature of 70° F, 
while moving with an angle of attack of 1 . 5 ° at a velocity of 253 feet 
per second through air at a uniform temperature of 3 O 0 F. The values of 
f Cx computed by the different methods for this example are shown 

graphically in figure 3 . 

It is seen that all of the methods discussed give answers that are in 
fairly good agreement. The reasons for differences are readily observed 
by inspection of table I in which the various terms of equations ( 23 )* ( 28 ) 
( 33 ) j and - (39) are presented. For example Allen and Look*s values are 
obviously too low because the calculations were based on a Pr value of 1 . 
The values calculated by equations ( 33 ) and (39) are in very good agreement 
but some of this agreement is fortuitous* as table I reveals. Thus* the 

8p 

individual corrections of Squire for — and Pr are of such a magnitude 

dx 

that their combined effect yields results which are in very close agreement 
with those of Martinelli and others; calculations on a different airfoil* 
however* may give results which are considerably more at variance. 

If, as Is often the case In practice, an approximate value of the 
external conductance* or a value indicating an order of magnitude* Is 
desired* then the method of Martinelli is the simplest and most rapid and 
Is usually of sufficient accuracy; In other cases the accuracy of the 
answer desired would determine the choice of method. It should again be 
remarked in this connection, however, that the example cited for purposes 
of comparison does not represent an extreme case of pressure gradient 
around an airfoil and that the deviations among the methods due to the 
effect of the pressure gradient on the heat transfer might be considerably 
greater for other types of shapes* such as thick wings, fuselages, and 
so forth. 


Methods Employed In the Turbulent Regime 

In order to compare the various methods for calculating the unity 
conductance along an airfoil for the turbulent boundary layer, the wing 
profile utilized by Frick and McCullough (NACA 65*2-016) was selected for 
calculation. The airfoil has a chord of 7 feet and was assumed to be 
moving with a velocity of 206 feet per second. The average temperature 
of the air in the ' turbulent boundary layer was taken as h0° F. 

The five methods of calculation previously discussed were applied to 
the airfoil, assuming h turbulent boundary layer to exist from x = 0. The 
results are plotted in figure h. It is noted that the five methods yield 
results, which, as in the case of the laminar— boundary— layer calculations. 
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are in fair agreement. Certain differences are apparent, however, which 
merit further discussion. 

• In order to facilitate comparison of the various methods, the pertinent 
tewns of equations (70), (73)> (7l)*(75)j and. (8l) are presented in table H. 
It is apparent from this table that Frick and McCullough's original 
equation is too low because of the tacit supposition that Pr = 1. Inclusion 
of the . approximate correction for Pr in Frick and McCullough's methods 
raises the curve about 20 percent. The remaining methods check closely 

up to — ■= 0.4. At this point Squire's heat-^balance method diverges 

c • _ 

■rapidly from the other curves. Experimental, evidence is not available to 
check this phenomenon, but it is probable that the heat-balance method 

proposed by Squire overemphasizes the importance of 
. > 

The lather close comparison, of the method proposed by Mart ine 111 and 
others with the more refined techniques of the other authors is partly 
fortuitous, since for airfoils with abrupt pressure gradients the results 
from the various methods may be considerably more at variance. 


Unit Conductance at the Stagnation Point 

Allen and Look and also Frick and McCullough suggest the use of the 
approximatl on 


6 


2 

stag 



5Re 0 c 


( 82 ) 


for computing the value of the point conductance at the_ stagnation point, 
where r Is the radius of curvature of the leading edge. This equation 
considers the airfoil leading edge to be elliptical in form. The equation 
of the point conductance at the stagnation point, then becomes 


Nu Btag = M = 2.42^ = 2.42/Ii^ (8 3 ) 

where d is twice the radius of curvature. 

At the stagnation point the equations of Squire reduce to the fom 

M 

&i 2 = ~ (84) 

Re D 
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Since 


f _ 0.5715 

c x " d 



(36) 


Hu stag = ~ = 1.63 \jBej) ( 


8, 


(85) 


8 1 

where — is found "by means of the equation 
8 2 


'Sri 2 J20 = 

, 8 1 / \ 8 i / Pr 


( 86 ) 


The tern — accounts "both for variations in Pr frcan. unity and for the 
80 

pressure gradient about the stagnation point. At Pr = 0.7 20 equ ation ( 86 ) 
8 n 

gives — * = O. 58 I so that 


^ u stag _ 0.95 


(87) 


The equation proposed hy Martins Hi for the stagnation-point conductance 
is that derived in the theoretical analysis of heat transfer from, the 
forward end of a cylinder, namely 


Nu sta& - l.lW^ep 0 ^ (29) 

or at Pr = 0.720 

^stag = 1-0 ^©d • (88) 

Comparison of equations ( 83 ), ( 87 )^ and ( 88 ) immediately indicates 
the degree of correspondence between each of the three methods and the 
stagnation-point conductance for the leading edge of a cylinder. In the 
method of Martinelli, the flow at the leading edge of the airfoil is 
initially postulated to be exactly that of the flow at the forward portion 
of a cylinder and hence the airfoil stagnation-point conductance corresponds 
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exactly with that of the cylinder. Apart from this method It is seen that 
Squire's method gives veiy much better agreement with data obtained at 
the stagnation point of a cylinder (reference 7, ch. IT) than does the 
method of Allen and Look. _ 

If the cylinder-stagnation-point value of f 0 is considered as being 

very close to the actual leading-edge value for the airfoil, then the error 

made by Allen and Look in setting the ratio -gi equal to unity is well 

“2 

illustrated by the wide discrepancy between their value of the point 
conductance at the stagnation point with that obtained from analysis of 
heat transfer at the stagnation point where the pressure gradient is 
extreme. This discrepancy resulting from use of the Allen and Look method 
is an indication of the superiority of Squire's method when applied over 
the entire airfoil. 

If in Allen and Look’s method the leading edge of the airfoil is 
considered to be cylindrical instead of elliptical, their equations reduce 
to 


2 0.289rc 

>Bta S " Be„ 


( 09 ) 


and 


®Wg - ¥ - 2 -o i/S» 


(90) 


rather than equations (82) and ( 83 ). 

When the expressions for the point conductance at the stagnation 
point are applied to the Joukowski profile, which was used as an example 
in comparing the methods for the laminar regime, there is obtained: 


Method 
Martins lli 

Allen and Look (elliptical) 
Allen and Look (cylindrical) 
Squire 


Cx stag 

(Btu/(hr)(sq ft)(°F)) 

109 

26 k 

218 


104 



NACA TIT No. 1453 


39 


Propeller Calculation 

The problem of heating propellers to prevent the formation of ice 
is becoming of increasing importance in the aircraft industry, and same 
method of calculating the point unit thermal conductance on the outer surface 
of the propeller would therefore be desirable. 

The flow conditions around a propeller are complex, and the theory is 
rather imperfectly developed so that an exact analysis of the problem in 
a manner similar to that' employed for an airfoil section in this report 
does not at present seem practicable. 

It is well known, however, that a propeller may be considered as 
being made up of a series of airfoil elements, and a rough approximation 
to the variation of the point conductance radially and chordwise might 
therefore be obtained by calculating the ohordwlse distribution by any of 
the methods described herein for an airfoil shape, for each of these elements. 
Such a procedure was adopted here by applying Martins 111 equation for the 
laminar and turbulent cases to four different sections of a propeller Bhape 
whose characteristics are described by Bairstow (reference 25, p. 664). 

Because experimental data on the pressure distribution around the 
blade elements of the propeller were lacking, the approximate equation of 
Seibert (reference 26) 


U = u (l± — 2l» ^ 

80 \ 4 cos a/ 

+ for "upper" surface of airfoil 
— for "lower" surface of airfoil 


(91) 


was used to calculate the velocity at the edge of the boundary layer. 

The propeller in question was assumed to have a radius of 5 feet, to 
be rotating at 2000 rpm, and to be maintained at a surface temperature of 
70 c F, in air at 30° F. The lift coefficients for each of the sections 
were obtained from reference 25 (p. 664) in which tables of data are 
reproduced which were taken at the National. Physical Laboratory. It was 
assumed that the angle of incidence of each of the blade elements was 
the same and equal to 6°. A plot of the results and a diagram of the 
propeller selected for the example are given in figures 5 sn d 6, respectively. 

In calculating the heat loss from the surface of the propeller, two 
cases were postulated: (l) The flow remains laminar for the entire chord 

of the airfoil section for both the upper and lower surfaces, and (2) the 
flow is turbulent over the entire chord for the upper surface and laminar 
for the lower surface. In each of these cases an average f c for each 

section was found by graphically integrating under the curves of figure 5 
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and the total heat loss was computed "by adding the heat dissipation 
from each of the four sections. In the turbulent case this total heat 
loss was found to be 20,487 Btu per hour, whereas in the laminar case a 
value of 5772 Btu per hour was obtained. 

Very little is known about the location of the transition point on 
propeller shapes, but it is probable that the flow is not completely 
laminar except possibly at positions close to the hub. Until further 
information is obtained regarding the nature of the flow about a propeller, 
it is probably safest to assume the turbulent value for the heat loss as a 
design estimate. 


Department of Mechanical Engineering 
University of California 

Berkeley, Calif., September 29, 1944 



APEEKDU A 


MOT® OF DEKCVAHON AHD PROCEDURE F03 CALCULAmSG POINT CONDUCTANCE 


The general energy equation for an ideal gas may he written (reference 7, vol. II, 

p. 607) 




(92) 


where D indicatee a total derivative and 4> is the "dissipation function." 
For the case of two-dimensional flow expression (92) reduces to the fom 


pjc„fS + u^ + T S'| _( , ^ + u Se + t 2£') + + . 

P H&t a* ay ) 1st c* Sy ] Va? ipj 


(93) 


Under the usual postulates of boundary-layer theory the terms ^ and are set equal 

a \ 2 SxC 

1 is substituted for $ 


to zero; when p. 


pJCp 


(gu ug + t » &V •* * (¥ 

ytt Sr bjJ \5t Si/ ay 2 \dy. 


‘) 


(94) 


If 6g is the thickness of the thermal boundary layer, on integrating equation (94) from 
y = 0 to y c 5]j, 


JC, 



„ dT . St &rV 

p — + pru — + pv — I dy - 
St Si byy 



— + u dy = -Jk 
St hxJ 


rs-r« 


dy (95) 


•fc" 

H 
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since 


and 


Also 




vanishes at 


y ® &B* 


which is the outer edge of the thermal layer. 


Now 


■p" 

ro 



pu 5 ^ = 5 



pvx * - r* * & <pu) 37 ■ (pui v-=% I s 

pv Saj.pvrl 511 - /** T g| (pv) dy 


(96) 


(97) 



vanishes for y = 0; and for y = 5^ it follows from the equation of continuity that 

(oT) y=% = ' [t - s (pu) ] 37 

1 i 

S pudy + (p)yD&H 5T + (pu, y=&H ST (99) 


Hence mating use of continuity there is obtained 


NACA TN No. 1453 





+ ~ P u3? dy _ Tqf^ p dy + pu dy ) (100) 


Sx 


Sx 


■ where is the temperature at the outside of the layer. When equation (100) is substituted in 

equation ( 95 ), the following equation for the integrated energy balance on the thermal boundary 
layer is obtained: < 


JC p 


St 



pT dy + 


S 



puT dy - Tq 






+ u 

Sx 



( 101 ) 


* * 


The second aforementioned general integral relation is the well-known von Karman 
"momentum equation," 


_S_ 

St 



pu dy + 


S_ 

Sx 



pu 2 dy 



p dy + 


S 

5 



5 


Sp 

£ 


-V 



\=0 


( 102 ) 


If steady flow and an incompressible fluid are postulated and the dissipation term is 
neglected, these two equations reduce respectively to 
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& l* * 37 - & j[ % u 37 - - a d) y=c (w3) 

sj[ uS ^- v ^f 0 udJ "-?s- v (sf ) y=0 < 104 > 


When a fourth-power polynomial Is postulated to represent the velocity 
distribution in the 'boundary layer, that is 

^ “ ay + by 2 + cy3 + dy^ ( 105 ) 

with the boundary conditions 

at y ** S: 


u = IT 


Su 


du 



U* 


0 

0 


at y = 0; 
u = 0 
v = 0 


therefore 

jfu = _1_ dp = _ UU» 
dy 2 pv dx v 

Pohlhausen’s expression (reference 27) is obtained 







( 106 ) 


( 107 ) 
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where 


X 

v dx 


( 108 ) 


Substitution of the drag at the wall from this polynomial into the 
momentum equation yields as an equation for X 

i* = ^ bM + [i 2 h(D + x] (109) 


where g(X) and h(X) are known tabulated functions of the argument 
(reference 'J, vol I, p. 160 ) , and primes denote differentiation with 

respect to x. At the point of stagnation X «= 7.052 and X* = 34 . 05^’: 

U 

with these initial conditions equation ( 109 ) can be integrated either 
numerically or by the method of isoclines. If the temperature distribution 
is assumed to be of the same general form as that of velocity, namely 



+ by 2 + 


oy3 + dy^ 


with the boundary conditions 


( 110 ) 


at y = 5^: 


T = Tp 


ST _ mt 

ax “ 


3T 

¥ 

ay 2 


0 

0 


at v = 0; 


u = 0 
v = 0 


therefore 


d%F 

§y 2 


0 


(in) 



there 1b obtained 


■t=* 

o\ 


I “ 2 d) " £ (4) 3 + (kf 

Equation (103) may he written in the following form, 

h OTi S £ * - T i s D j[ % ^ - - a d) 


( 112 ) 


(1130 


If the "kinetic temperature rise" is set equal to zero, then the temperature at the edge of 
the boundary layer ie the same as that in the free stream, that is = 1^, and 


A. 

ax 


'5g 

Dl I 




^ dy 
U ^ 


-£(SL 


(114) 


Substituting equations (107) and (112) and integrating, 




(1 + O.O83X) — + (-O.500 + 


( &rr\3 

0.1251)1-=/ + (0.200 - 0.0331) Ur J 
\ o / \ 0/ 


u JD 


ay 3 *h 


■1 i • . 

(0.866 + 0.0721)^ + (-0.479 + o.H9i)^) 3 + (0.194 - 0.0321)^ 


(115) 


(116) 


It is well known that as a good approximation (reference 7) 


^3 z Pr 
& 


(H7) 
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bo that we may define two functions 

J ^ g JL ay = F^l,!*) s (0.866 + 0.072l)Pr -1 / 3 

+ (-0.479 + 0.119X)Pr -1 + (0.194 - 0.032\)Pr"^/ 3 (118) 

V & ~ F 2 (4,Pr) h (1 + 0.083X)Pr -1 / 3 +- (-O.5OO + 0.125X)Pr“ 1 

+ (0.200 - 0.033 4)Pr"V 3 ( 119 ) 

and let 

Cp S3 Fj|^ Fg 

Then equation (113) reduces to 

(ieo) 

Since all temperatures are measured with the surface temperature as a 
datum, 

f c x = rc p ^ (n&gq)) x 3600 (121) 

Q3C 

The procedure for calculating the point conductance in ary given case 
would then be: 

(1) Calculate U' and ■C’* either graphically or analytically. . 

(2) Compute the ohordvise distribution of X by Pohlhausen’s method 
from the experimental pressure distribution. (See equation (57a).) 

(3) Find 6 by means of the relation X = ^ U*. 

(4) Find 8 jj by means of the relation = Pr —1 /3 # 

5 

(5) Compute the functions F^ and F 2 and hence cp. 

(6) Substitute in equation (121) and evaluate the derivative graphically. 



It is seen that this procedure permits extension to cases in which compressibility and 

dissipation effects are to "be considered. "When p = ^ and - -pU0 f are substituted 

xu OCC 

in the general energy equation for steady flew. 


JC 


p 5 


r* pui * ' Ti s r •>“ * ) -r u i * ■ * “ r $) a * ^ 


then 

Jsir*— -^.Ts.) * s r J t>- fdj 


(123) 


JC, 


ls ,p 


f- 


ay I -h Aip 

^ 1 E dx 1 * 



% 


J dy 

T 




(I2t) 


ii P 

R ^ Jo 


^ Jo 


R a* Jo T 


S521 


•fA 5»--*(fL 


•m* 


(125) 


£ 


K+ri ’Ofr hi vovm 



", r«a. _ » / r^a \ , „ r®H „ t, A /,r ®, \ 

^‘J, S* + 3s('J, ^J+S'™}, S**-TfE^J, «**/J 


pn* u T i , - /St 

+ rr -rr- dy “ — Jkf ■c— 

®U0 u T (*. 


r% 

w 


( 126 ) 


[**«*& r* s ** e £- Fn«*^r- 


\*‘ -1 


^ a a r«u^i ^ /■*/ py 

EH Jo O * >i»il ” V/pO Jo V 5 ?/ 

The polynomials of equations (107) and (112) permit the evaluation of the integrals 

P 5 H u r 8 a n T 1 [ S H~ U ,2 

- dv. I - dy. and I ( t dLv 

Jo U Jo U T Jo W ' 


“ -x J'- 


Pb . n 


(127) 


(128) 


When these integrals are called. F-^XjPr), FgfX^Pr), and F^X^Pr), respectively, equation (127) 


■becomes 


( J. + JJt\ 6JT, +( s5n - u* - 53 l 1 + Jl-^&nFo + uZ&trFi ‘ + &rr*Fi ^ 

E UETi ; u x V ETx T X JC p T / “ c \ “ x “ V 

.1 

- ^l^V 5 ^ + & H F 2 , )J ; " " Jk (^) yr<) + ^3 
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C p p U*(l + BTi) 
E ETp 

u(6h‘Fi + &#!*) 


BgFi + 


iAj* 

C™ i 


U*Ti* U' 
- U» i- + 


<JC,,T 


B#2 


uti 2 ( &hF 2 + BaFg* ) 


+ HI3 


(130) 
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APPENDIX B 

DISCUSSION OP CERTAIN BOUNDAILT— LAYER CONCEPTS 


In the development of boundary— layer theory, certain basic postulates 
are tacitly made which, because they are not clearly stated, may cause 
the engineer who is not primarily an aerodynamic is t considerable difficulty. 
The purpose of this appendix is to discuss several of these points somewhat 
more thoroughly than is done in the usual references on aerodynamics. 

There are two general methods of analyzing boundary— layer problems: 

(a) The solution of the boundary— layer differential equations 
(equation ( 132 )). 

(b) The solution of the boundary— layer momentum equation 
( equation ( 104 ) ) . 

The first part of the appendix deals with method (a) ; the second part 
deals with method (b ) . 


Boundary-Layer Differential Equations 

The incompressible— flow bound ary— layer differential equations are 
derived from the general hydrodynamic equations with the following 
postulates : 

(a) The flow pattern about the object is two— dimensional. 

(b) The thickness of the region next to the solid surface, in which 
the velocity gradient is large, is small compared with the other linear 
dimensions of the object. 

(c) The flow is incompressible. 

(d) There is no separation of the flow from, the solid object. 

(e) The fluid in contact with the solid surface has no velocity 
relative to that surface. 

Reference 7 (vol. II, p. 610) presents a derivation of the boundary- 
layer equations, for curved sxirfaces, on the basis of these postulates. 
The boundary— layer equations are: (See also fig. 1.) 



52 


NACA TN No. 1453 



^ + ^1 = 0 
5x by 


(131) 


The following further postulates are visually made in solving 
equation ( 131 ): 


(e) Telocity u is not a function of time. 



Is negligibly small even along a curved surface, so 


that the pressure is postulated to be invariable in the direction normal 
to the surface. 

The equations then reduce to: 


u 


du 

dx 


+ v 


du 

by 


1 5p 
p dx 





du | bv 
dx by 


m 0 


(132) 


By making postulate (f ), all problems of two-dimensional boundary- 
layer theory beocme problems of flow along a plane surface (flat plate) 

along which the pressure gradient : varies in seme given manner. 

The boundary conditions imposed on the bo undar y— layer equations, 

3p 

regardless of the variation of ^ 


are 
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at y = -0 
and x a oo 



0 

0 


at y 


u = U 


and x = 0 


M » o 
°y 

d 2 ^ 0 

5? 0 


(133) 


The extension of the solution to y = « assumes that the whole fluid 
field is in viscous motion. With these boundary conditions, the velocity 
profile at any x will have the form shown in figure 8. 


Since at y 


00 , 


U = U, 



equations far from the plate 1)6000161 



the boundary— layer 


i sp tt au 

- D S 


(13M 


Integration of equation (13*0 yields Bernoulli. *s equation. 

Since ^2 does not vary with y, the boundary— layer equations 
become: 


u 


cm 

3x 


+ v 


5u 

Sy 


= U 


8v 


8u 

3x 


+ v 






= o 


(135) 


As was previously mentioned, these equations are strictly applicable 
to a flat plate only. Approximate solutions of these equations over curved 
shapes are obtained by postulating various variations of XJ with x and 

r¥TT 

'substituting this resulting form of U ~ into the foregoing equations. 

“ax 

Only certain special forms of U = f( z ) can be handled mathematically . 


The following table was compiled frcm the solutions of the boundary- 
layer equations presented in reference J r and reveals how flow about 
various objects is identified with certain variations of IT along a plane 
surface. 
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u 

_ i Sp D y au 
p cix 3x 

"Equivalent flow system" 

Constant 

0 

Flow along flat plate 

P x x 

Pi 2 * 

Flow near stagnation 
point of cylinder 

c. 

X 

c 2 

x3 

Flow in converging 
channel 

' ex ' 11 

c^m x 2m_1 

m = 0 Flat plate 
m = 1 Stagnation point 
m = — 1 Converging channel 

m=oo 

51 ^ 

m=l 

1 1 

Dl=oo 

5 o 2 m x 2 ^^- 
m=l 

Any curved shape if series 
is known from experimental 
data 


A diagram of the velocity distribution near the stagnation point of 
a cylinder, obtained from a solution of the boundary— layer equations, is 
shown in figure 9 . A discrepancy with the actual physical system is 
immediately apparent. Figure 9 shows the velocity far from the object 
increasing linearly with x, whereas obviously in the physical system 
the velocity far from the surface remains constant. 

The flow pattern about the actual physical system is shown in 
figure 10. This figure reveals that the velocity in the physical system, 
at any fixed value of x, increases as y increases, reaches a maximum, 
and then decreases to an asymptotic magnitude as y -» ». Comparison 

of the velocity distribution about a cylinder obtained from potential 
theory (zero viscosity) with the experimental points would show a good 
check until the surface of the cylinder was approached closely; the 
potential solutions would then became greatly in error. Conversely the 
solution of the boundary— layer equations would show good agreement with 
the data for magnitudes of y from zero to about- the point of maximum 
velocity; beyond this point the boundary— layer solution would deviate 
greatly frqm the data. 

Because of the deviation of the boundary— layer solution from experimental 
results for flow around curved objects, a limit crust be placed on the region 
in which the boundary— layer solution is applicable. This limit is called 
the thickness of the boundary layer. In analyzing boundary-layer problems 
aerodynamicists as a rule concern themselves with what occurs within this 
boun dar y layer and usually neglect completely the flow outside of this 
boundary layer. 
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It is -evident that a unique definition of the thickness of the 
■boundary layer o is difficult to' establish. Several possible definitions 
are: 

(a) The point y, where the boundary layer and potential solutions 
intersect 

(b) The point y, where u = O. 99 U, or any other arbitrary fraction 
of TT 

(c) The point y at which the total pressure reaches a fraction F 
of the total pressure in the free stream. 


Definition (c) of S (for incompressible flow) is equivalent to the point y 
at which u = /F U. 


(d) A characteristic length, called the displacement thickness 
Sj_ may be defined, so that 


U dy - (U&i) 



(136) 


If the flow along the body obeyed the boundary— layer equation, the rate 


of flow of fluid an rose any x 


would bo 



This is less than the 


n 00 

quantity / U dy because of the retardation of the flow near the surface 

Jo 

of the object. The difference is called TJo^, thus defining the character- 
istic length 8 j_. The length &i is not the boundary— layer thickness, in 
itself, but is related to it. 


The exact definition of boundary— layer thickness adopted should not 
influence the final results of a boundary— layer analysis as long as the 
definition is utilized in a consistent manner, since the analytical so lu tions 
of the boundary— layer equations never involve the boundary— layer thic kne ss 
directly. 


From figure 10 it is evident that in the physical system the velocity 
outside the boundary layer varies with y. In aerodynamic analysis, 
however, the velocity U Is called the velocity at the edge of the 
boundary layer. The velocity XJ is defined by 


dtf _ _ 1 dp 
dx p dx 


(3.34a.) 
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whence 


— = - - p + Constant 
2 P 


P + 



Constant 


Thus U may be calculated from the pressure distribution around a curved 
object. It is defined by equation (134a) and is the velocity just outside 
the bou ndar y layer only in the sense that for any fixed value of x, at 
one point y outside the boundary layer, the velocity U will exist. The 
correct value for the velocity at the edge of the boundary layer will, of 
course, depend upon tho definition of the boundary layer adopted. Utilizing 
U indiscriminately for the velocity at the edge of the boundary layer in 
conjunction with different definitions of the boundary— layer thickness will 
lead to erroneous results. 

It has been shown that the boundary— layer* solution 1 b not applicable 
outside the boundary layer. If now it is postulated that the flow outside 
the boundary layer does not affect the velocity inside the boundary layer, 
the laminar-boundary— layer equations may be applied to a Bystem even 
though the flow outside the boundary layer is turbulent, as long as the 
flow within the boundary layer remains lam inar. Once the flow within the 
boundary layer becomes even partly turbulent, the boundary— layer equations 
are no longer applicable. 


Momentum Equations 

By the considerations of the flow within the boundary layer it may 
be established that a certain thickness 5 exists in which the boundary- 
layer equations are applicable. If it is furtlier postulated that: 

(a) the exact form of the velocity distribution outside the boundary layer 
Is of no Importance whatsoever in determining the behavior inside the 

boundary layer and (b) that a velocity U calculated from p + p = Constant 

exists at the edge of the boundary layer, certain simplified equations may 
be written which allow the approximate analysis of the velocity distribution 
within the boundary layer. 

Since the two postulates mentioned are not exactly true, the 
momentum equations are only approximations . Experimental results indicate 
that tho approximation is fairly good. The momentum equation for steady 
flow (see appendix A) is 


_b 

dx 


( 104 ) 
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It should, he particularly noted, that integrations extend, only to 8 instead, 
of co as was the oase for the houndary— layer equations. Thus the "boundary- 
layer thickness is- of primary importance in the momentum equations, although 
its definition was really not essential in the analytical solution of the 
houndary— layer equations. 

The momentum equations are solved hy expressing the velocity ratio ^ 
as a function of Z, A typical expression for H (appendix A) is 

5 - <!) + Kl) 2 + °(l) 3 + Kl) 4 


Since only the region between — = 0 and — = 1 is being considered, 
the following approximate boundary conditions may he imposed on — 


0 

II 

at 

1 = 1 

y 


8 

0 

n 

<1 

at 

t m 1 


8 

0 

11 

*1 

at 

1 = 1 

v-3 

8 


u 

at 

£ = 

TJ 

8 

5u 

8y 

is finite at 

y _ 
8 

dfu 

= 5 at 

z = 

Sy 2 

X 

8 

S3u 

= 0 .at 

Z = 

3y3 


8 


(137) 


(See reference 7, vol. I, p. 156 ). These boundary conditions, which are 
not exactly true in the physical system, allow, however, the evaluation 
of the constants a, h, c, d, and so forth. Obviously the velocity dis- 
tribution is only applicable between ^ = 0 and ^ = 1. Outside ^ = 1 
there is no relation between the polynomial expressing ^ and the 
experimental velocity distribution. 

U. Y 

Having ^ as a function of the momentum equation may be solved 

for the shear at the wall, the boundary— layer thickness 5, and the 
displacement thickness 8 ^. In addition, for purposes of analysis the 
momentum thickness defined by 

PS 

9V 2 = (U - u)u dy (138) 

Jo 


is often calculated by the momentum, equation method 
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The advantages of the mcBientum method are the relative simplicity of 
accounting for variable pressure gradients as a function of x and the 
fact that it may "be applied approximately to both laminar and turbulent 
boundary layers. 
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APPENDIX C 

SOME REMARKS ON TEE SIMILARITY OF VELOCITY AND TEMPERATURE 
DISTRIBUTIONS IN THE PRESENCE OF A PRESSURE GRADIENT 


As was discussed in the section of the text LAMINAR REGIME, Method 
of Squire, inspection of the boundary— layer equations, for both lami nar 
and turbulent boundary layers, reveals that within the boundary layer the 
pressure gradient will probably cause only secondary differences in the 
temperature distribution as compared with the velocity distribution. 

This conclusion is even more important for the region outside the. boundary- 
layer. As shown In figure 7 and as discussed ih appendix B, the solution 
of the boundary— layer equations, even in the presence of a pressure 
gradient, yields values of U and T which asymptotically approach 
certain magnitudes of U and T which are supposed to exist far from, 
the solid— fluid interface. 

Physically, however, in the presence of a pressure gradient the 
actual velocity outside of the boundary layer does not approach the velocity 
U asymptotically, but (in a region of negative pressure gradient) reaches 
a maximum and then decreases to a magnitude of the free— stream velocity u«j. 
The temperature distribution, however, (neglecting frictional heating) 
approaches T M asymptotically. 

It Is apparent therefore that, at least outside the boundary layer, 
the pressure gradient affects the velocity distribution much more than 
it affects the temperature distribution. 
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TABJE I 

CALCULATIONS FOR LAMINAR BOUNDARY LAYER FOR JOUEOWSKI_ AIRFOIL 
[o = 7-78 ft, u*, = 253 ft/sec, 7 =- O.O 78 lb/ft-3, 

Cp = 0.24 Btu/ (lb) (°F), Pr = 0.?2, a = I. 5 0 ] 


x/c 

IT 

u 

•0 

Re x 

0.332 
V B 0 x 

[A ]' 1 / 2 

( 1 ) 

( 2 ) 

Pr-2/3 


0.06 

1.270 

76.4 x 10 ^ 

3.80 x io - ^ 

1.197 

1.059 : 

1.242 

1.22 

.10 

1.276 

131.5 


2.81 

1.121 

1.050 

1'. 242 

1.24 

.16 

1.276 

209 


2.22 

1.Q70 

1.012 

1.242 

1.24 

.20 

1.272 

263 


1.99 

1.042 

1.012 

1.242 

1.24 

.28 

1.264 

366 


1.68 

1.003 

•997 

1.242 

1.25 

•36 

I .251 

466 


1.49 

.968 

.980 

1.242 

1.26 

. 44 

1.229 

559 


1-37 

'.910 

•959 - 

1.242 

1.28 

• 52 

1.218 

655 


1.26 

.895 

• 932 - 

1.242 

1.30 

.60 

1.198 

744 


1.18 

.845 

.895 - 

1.242 

1.32 

;72 

1.162 

860 


l.l4 

• 770 

.848 

1.242 

1.36 





f 0 as .determined by - 

- 




Allen and Look 

Frick and 

Martinelli „ 

Squire 

x/c 

(equatkan ( 23 )) 

McCullough 

and others 

(equation (39)) 


(3) 

(equation ( 28 )) 

(equation (33)) 



0.06 

9.8 



12.2 

10 

.2 

10.6 

.10 

6.8 



8.50 

7-59 

7-95 

.16 

5.16 


6.42 

5.98 

6. 

05 

.20 

4.51 


5-6o 

5.38 

5-43 

.28 

3.65 


4.54 

4.51 

4.54 

•36 

3.08 


3.83 

3-94 

3.96 

.44 

2.60 


3-24 

3 

57 

3-51 

•52 

2-35 


2.94 

3 

27 

3.18 

.60 

2.04 


2-53 

3-00 

2. 

86 

.72 

1.74 


2.16 

2.82 

2 - 61 ...1 



3 


Calculations based on. 


Pr value of 1 






wnn n 


cAicramaira foe otrbolhti bockdabi iaue for haca 65,2-016 airfoil 

£ calculated In reference llj + lje approximately calculated by utilizing | 

a. 0.0296 EBj 0,2 B'/5 - ^ 

*- ^3 7^ *WT 


0.0296T 


14.0 x 10^ 

40.2 

65.2 

11.3 

13.6. 

25.2 

36,5 

47.9 

50.4 


Martinelll 
and. cHierei 
(equation (71)) 


0,00277 

.00224 

.00204 

. 00102 . 

.00174 

.OOI55 

,ooi44 

.00136 

.00131 


-2.2 x io _q 

-3-74 

-3-54 

-2.94 

-1.80. 

-1-71 


a a determined by - 


Frlok end 
MftOnl 1 micft 

I 

(aquation (70)) 

Frlok and 
MoCiiU miph 

n 

(equation (73)) 

Sqjllre 

I 

(equation (75)) 

Squire 

IX 

(equation (60.) ) 

68.5 

84.4 

80.4 

73*4 

42.4 

52.6 

48.4 

50.6 

38.5 

47.8 

43.6 

46.0 

33.2 

41.3 

37.7 

39.2 

32.4 

40.6 

36.4 

36.0 

28.8 

35-6 

32.3 

32.2 

26.7 

33.3 

29-9 

20.8 

25.4 

33-5 

31.6 

29.1 

28.4 

26.3 

30.8 


5 J(Pr - 1) + log, [l + 5 fPr - l)]! 
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Figure 2.- Joiikowski profile used in comparison of methods for 

laminar regime. 


5/a = (cos e + 0.04) 


(cos 9 + 0.04 ) 2 + (sine + 0.05)^ 


n/a = (sin 9 + 0.05) 



083 ; 

(cose + 0.04 ) 2 + (sin e + 0.05) 2 
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Figure 3.- Comparison of methods for calculating point unit thermal 

conductance in laminar regime for an airfoil section (Joukowski profile). 
True airspeed, 253 feet per second; air temperature, 30° F; wing 
temperature, 70° F; angle of attack, 1.5°; wing chord, 7.78 feet. 
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x/c 

Figure 4.- Comparison of methods for calculating point unit thermal 
conductance in turbulent regime for airfoil section (NACA 65,2-016). 
True airspeed, 206 feet per second; air temperature, 30° F; 
wing temperature, 70° F; lift coefficient, 0.55; wing chord, 7.0 feet. 


rve Section 

l F 

3 E 

: d 

) c 

(See fig. 6) 


Laminar 


Turbulent 






x/c - Lower surface 


x/c - Upper surface 


Figure 5.- Point unit thermal conductances on lower and upper surfaces of a propeller. Radius, 

5 feet; speed, 2000 rpm; angle of attack, 6°; surface temperature, 70 0 F; air temperature, 30° F. 
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r = 0.412 R, Chord = 0.164 R, C L = 1.172 



r = 0.602 R, Chord = 0.602 R, C L = 1.141 



r = 0.75 R, Chord = 0.137 R, C L = 1.061 



r = 0.88 R, Chord = 0.137 R, C L = 0.961 
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pressure gradient. 








